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Abstract: 

The fractions of samples spanning a lattice at its percolation threshold are 
found by computer simulation of random site-percolation in two- and three- 
dimensional hypercubic lattices using different boundary conditions. As a 
byproduct we find Pc = 0.311605(5) in the cubic lattice. 
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Introduction: 

In the case of random site percolation, the probabihty that the system is 
percolating below the percolation threshold {pc) is zero and it is unity above 
Pc- What is the probability of a cluster to span the lattice at the percolation 
threshold [Q? The earlier studies [2-6] tried to establish this quantity (span- 
ning probability at percolation threshold) as an universal quantity. To check 
the robustness of the spanning probability, we have studied here the effect 
of different boundary conditions on the spanning probability at percolation 
threshold. Our study is confined to two and three dimensions and restricted 
to random site percolation. 

For free boundary conditions, in three dimensions, the critical spanning 
probability is 0.265 ± 0.005 and it is universal |^. However, an earlier study 
[0] gives different result (~ 0.42). This is due to the different boundary 
conditions used by two groups of researchers as pointed out in [4,6]. Here, 
we study the critical spanning probability (in two and three dimensions for 
random site percolation) and its dependence on the different boundary con- 
dition. 

For clarity, let us first define our boundary conditions: In D-dimensions 
the spanning direction is always free 

(1) Free boundary condition (FBC): The other D-1 directions are also free. 
(2) Helical boundary condition (HBC): The Helical boundary conditions are 
used in all D-1 remaining directions. (3) Mixed boundary condition (MBC): 
This is not applicable in 2 dimensions. In 3 dimensions, the Helical boundary 
condition is only applied in one of the directions of the planes and the other 
direction is free. 

We are considering only site percolation. The microscopic connectivity 
is considered by the nearest neighbour occupied sites only and the macro- 
scopic connectivity (spanning) is checked by standard Hoshen-Kopelman [|I| 
algorithm. The spanning is checked only in one direction. 

Results: 

In D = 2 we have taken p^ = 0.592746 For a fixed system size L we 
occupied the lattice randomly with probability Pc- We have used Kirkpatrick- 
StoU R250 random number generator as well as multiplication with 16807 [|^]. 
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Then we have checked whether the lattice is percolating or not by standard 
Hoshen-Kopelman algorithm for Ng number of different random samples. 
The number (fraction) of samples percolating at Pc is called the critical span- 
ning probability {Rp). For a particular boundary condition, the spanning 
probability Rp at Pc has been calculated for different system sizes (for exam- 
ple, ranging from L = 17 to L = 6551 in two dimensions). We have calculated 
the critical spanning probability in 2 and 3 dimensions for the maximum al- 
lowable (in 128 MByte CRAY) system size. The spanning probability, Rp 
(in two dimensions) has been plotted against L for free boundary conditions 
(see Fig. 1). From the extrapolation, we found Rp = 0.501 ± 0.003. This 
is in good agreement with the earlier result {Rp = 1/2) The spanning 
probability for helical boundary condition in two dimensions was given ear- 
lier as 0.637 and we report Rp ~ 0.638 in the asymptotic limit (L — oo) 
(Fig.2). 

In D = 3, we have calculated the spanning probability i?p at = 0.3116 
[H]. Here, also we have calculated Rp for different boundary conditions in the 
same way described above. Fig. 3 shows the plots of Rp against L for free 
boundary conditions. From extrapolation we obtained, Rp ~ 0.28, which is 
slightly higher than recently reported value. Using the mixed boundary 
condition, we have calculated the critical spanning probability Rp ~ 0.41 in 
the asymptotic limit (Fig. 4). This result is slightly smaller than the earlier 
estimate In the earlier case the mixed boundary condition is used 
in the following way: the top layer is fully occupied, the helical boundary 
condition is used in one direction of the plane. But in the present case the 
top layer is randomly occupied with probability Pc, and the helical bounadry 
condition is used in one direction of the plane. These two slightly different 
boundary conditions give different result in the smaller systems, however in 
the L — i> oo limit both boundary conditions give same result ^ . In addition, 
in three dimensions, we have calculated the critical spanning probability 
by using helical boundary condition in both directions of the planes. Our 
estimated value of Rp using helical boundary condition is approximately 0.513 
(Fig. 5). In each case, in three dimensions, we observed that the value of Rp 
is quite sensitive near pc- For this reason we have shown some results by 
taking pc = 0.31161. All these results also indicate that Pc lies in between 
0.3116 and 0.31161, compatible with earlier estimates p. 

In the present study, most of the data was obtained by taking A^^ = 10000 
in each processor of a CRAY-T3E computer (using 64 processors) . The CPU 
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time required, for the whole study, can be calculated from the following ex- 
ample: for system size L = 481 in three dimensions, each processor averaged 
over Ng = 800 different random realisations. The total CPU time required 
for this is approximately 11 hours per processor. 

Summary: 

In summary, we studied the critical spanning probability in two- and 
three- dimensional system using different boundary conditions. Our results 
from numerical simulations (with larger system sizes and using better statis- 
tics) show that the critical spanning probability is strongly dependent on 
the boundary condition. There are different alternatives to helical boundary 
conditions in 3 dimensions, as pointed out in 0. The value of the criti- 
cal spanning probability is quite sensitive to the value of Pc used. This has 
been shown in three dimensions and the size dependence of critical spanning 
probability can be an alternative tool for estimating pc- 
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Fig. 1. Size dependence of critical spanning probability (using = 0.592746) 
in two dimensions for free boundary conditions. Dotted line is a guide to the 
eye. 
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Fig. 2. As Fig. 1, for helical boundary conditions. 
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Fig. 3. Size dependence of critical spanning probability in three dimensions 
for free boundary conditions. The symbols • represent Pc = 0.3116 and A 
represent Pc = 0.31161. Dotted lines are guides to the eye. 
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Fig. 4. As Fig. 3, for mixed boundary conditions. 
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Fig. 5. As Fig. 3, for lielical boundary conditions. 
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